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THE FLUX-DEPENDENT RMS VARIABILITY OF X-RAY BINARIES IN THE OPTICAL 
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ABSTRACT 

A linear relation between absolute rms variability and flux in X-ray observations of compact accreting sources 
has recently been identified. Such a relation suggests that X-ray lightcurves are non-linear and composed of a 
lognormal distribution of fluxes. Here, a first investigation of the optical rms vs. flux behavior in X-ray binaries is 
presented. Fast timing data on three binaries in the X-ray low/hard state are examined. These are XTE Jl 1 18H-480, 
GX 339-4 and SWIFT J1753. 5-0129 - all show aperiodic (non-reprocessed) optical fluctuation components. 
Optical rms amplitude is found to increase with flux in all sources. A linear fit results in a positive offset along the 
flux axis, for most frequency ranges investigated. The X-ray and optical relation slopes track the source fractional 
variability amplitudes. This is especially clear in the case of GX 339^, which has the largest optical variance of 
the three targets. Non-linearity is supported in all cases by the fact that flux distributions of the optical lightcurves 
are better described with a lognormal function than a simple gaussian. Significant scatter around linearity is found 
in the relation for the two sources with lower optical variability amplitude, though observational biases may well 
contribute to this. Implications for accretion models are discussed, and the need for long well-sampled optical 
lightcurves is emphasized. 

Subject headings: accretion: stars - individual: XTE Jl 1 18H-480, GX 339^, SWIFT J1753.5-0129 - stars: 
X-rays: binaries - stars: optical: variable - black holes 



1. INTRODUCTION 

Several recent works have shown that the X-ray timing char- 
acteristics of accreting sources including Galactic black hole 
candidates, neutron stars as well as active galactic nuclei fol- 
low a linear relationship between their root-mean-square (rms) 
variation and average flux over a wide range of timescales 
(e.g. [Ut tley & McHardv 2001 [hereafter, UOl], Uttley 2004, 
iGleissner et al.l200 4 [hereafter, G04]). Such a relation suggests 
non-stationarity of the lightcurves; specifically, the level of in- 
stantaneous variance in a lightcurve i s not constan t , but i s rather 
linked to longer-term flux averages. lUttlev et al.1 (l2005l) [here- 
after, U05] point out that these properties imply a lognormal 
distribution of instantaneous flare strength. This, in turn, can 
be explained if the variations are composed not by a superpo- 
sition of independent shots, but instead are the result of accre- 
tion rate fluctuations operating at every location in the accreting 
flow of matter These fluctuations propagate inwards through 
the flow, effectively coupling perturbations at all inner radii. 
Interestingly, a similar 'rms-flux relation' has been shown to 
apply in other sources as we ll, for instance in a narrow-line 
Seyfert galaxy (Gaskell 2004) and in solar coronal X-ray and 
radio flares (Zhang 2007; Wang e t aL 2008). The relation thus 
seems to provide strong constraints on the fluctuating condi- 
tions of hot plasmas universally. 

In this Letter, a first investigation of the rms vs. flux behav- 
ior of X-ray binaries (XRBs) in the optical is presented. The 
presence of fast, aperiodic optical activity in so me accreting 
sourc es has been known since very early days (e.g. iMotch et aU 
Il982h . Recent work has found intriguing correlations be- 
tween the X-ray and optical variations in at leas t three X-ray 
binaries: XTE J111 8-H480 (iKanbach et al.H 200li) GX 339-4 
(lGandhietal. l l2008h and SWIFT J1753.5-0129 (iDurant et al.1 
120081) . The main driver of the rapid optical variability remains 
unclear, but several pieces of evidence disfavor an origin as re- 
processing of, or irradiation by primary higher-energy photons. 
Instead, it is likely that both the optical and X-ray emitting com- 



ponents share some common underlying cause, with exchange 
of energy between the components creat ing the complex in - 
terplay of timing patterns observed (cf. iMalzac et al.l l2004b . 
Studying the optical behavior in the rms-flux plane should pro- 
vide new clues on the underlying source of variability, inde- 
pendent of commonly-used correlation function analyses. The 
only other investigation of such optical behavior (as far as the 
author is aware), is for the a ctive galaxy NGC 4151 for which 
iLvutvi & Oknvanskiil (Il987h discovered a linear rms-flux rela- 
tionship in the U band. 

2. OBJECTS AND DATA 

The sources chosen for this study include all three XRBs 
known to show direct, accretion-driven rapid variability in the 
optical. The optic al data com e from the Skinakas/OPTIMA 
photometer (,Straubmeier et al.1.2001.) f or XTE Jl 1 18+480 and 
from the VLT/ULTRACAM camera dDhillon et al.ll2007h for 
GX 339-4 and SWIFT J1753.5-0129. They were observed at 
optical flux levels several magnitudes brighter than in quies- 
cence so the companion stars contribute negligibly in all cases. 
The most relevant details including used filter, time resolution 
and observation period are listed in Table [1] 

The corresponding X-ray data are from the RXTE/PCA de- 
tector, observed simultaneously with the optical (but see be- 
low). The RXTE ObsID prefixes are 50407 for XTE Jill 8H-480 
and 93119 for both the remaining sources, respectively. The 
counts and background were extracted over 2^4 keV for 
XTE J1118H-480 and over the fufl PCA energy range for the 
remaining targets, but the results herein are relatively insensi- 
tive to the exact range. Fuller details regarding the X-ray data 
can be found in the relevant papers cited in the previous section. 

The optical and X-ray data used for XTE J1118H-480 are 
strictly-simultaneous. For GX 339-4, RXTE monitoring ex- 
tends a few minutes longer on either side of the optical moni- 
toring; this full 62 min long dataset is used for better statistics. 
For SWIFT J1753.5-0129, changeable weather resulted in the 
best datasets being non-simultaneous; the X-ray data are 54 min 
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long and were observed on 2007 Jun 14 UT, while the optical 
data are from four days later (see Table [T]i. All sources main- 
tained a low/hard state throughout the optical/X-ray observing 
periods. 

Multiple nights of observation are available for each source. 
But a single night stood out as having the best optical weather 
conditions (in terms of wind for XTE Jl 118+480 and cloud 
cover for the other two) and only these are used herein. It is 
noted for completeness that using the strictly-simultaneous op- 
tical and X-ray datasets with worse weather conditions, i.e. the 
non-photometric nights of SWIFT J1753.5-0129 and GX 339- 
4, results in similar behavior in terms of the rms-flux relation 
described below, but with reduced significance. On the other 
hand, the nights of XTE Jl 118+480 affected by wind show a 
departure from this behavior; it is speculated that telescope vi- 
brations may cause loss of flux from the optical fibers of the 
OPTIMA instrument, rendering these data less useful. 

Power spectral densities (PSDs) for the sources in the liste d 
order in Table [T] can be found in ISpruit & KanbachI (|2002|) . 
Gandhi et al. (in prep.) and iDurant et alj ( 20091) . respec- 
tively. These have been recomputed here from averages of long 
lightcurve segments (256 s) of the same datasets used for the 
rms-flux relation below, and are displayed in Fig.[T] The stan- 
dard rms^-normalization has been used, so that integration over 

S iositive frequencies results in the fractional lightcurve variance 
Belloni&Hasingeij|l990h . 

3. COMPUTING RMS VS. FLUX 

To compute the absolute rms (cr) as a function of flux, I fol- 
lowed the procedure described by G04 (which is largely simi- 
lar to that described by UOl). In short, the lightcurve for each 
source was split into many segments of equal length ranging 
over 1-16 s, depending on the frequency range of interest (see 
below). The mean source flux in each segment was computed 
and this distribution was divided into a number of equal-size 
flux bins (F). The binning was optimally chosen so that the 
fraction of bins containing more than an adopted threshold of 
20 segments was maximized; only bins above this threshold 
were retained in the end. The PSD of each segment (rms'^- 
normalized) was computed, as was the expected white noise 
level. These two quantities were separately averaged, giving 
an average PSD in each flux bin as well as an average noise 
level. Finally, the mean level of the PSD over any particular 
range of frequencies (Az/) was computed, and white noise sub- 
tracted from this. This is the excess (above Poisson) mean frac- 
tional variance {var) in each flux bin, per unit frequency, cr 
is then simply \^vafAi' x F. Uncertainties on the mean a 
value in each bin are computed with usual periodogram statis- 
tics jvan der Klis 1989) and error propagation. 

Amongst the observations analyzed herein, there exists a 
wide range of source count rates, lightcurve sampling (AT 
and T values in Table [Til and intrinsic source variability lev- 
els (Fig.[TJ. So there is only a restricted selection of frequency 
ranges which are ideally suited for inter-comparison of the rms- 
flux relation. A base range of 0.5-5 Hz was found to be a good 
compromise for all sources in both X-rays and optical, also 
helping us to probe the broad-band optical noise continuum free 
from sharp features (e.g. near -7-10 Hz) for XTE Jl 1 18+480 
and from the low-frequency peaks below 0.1 Hz for the other 
sources; Fig. [T]). The rms-flux behavior over this frequency 
range is plotted for all sources in Fig. |2l For more complete 
comparison, the behavior over two more ranges is also plotted 



where possible: an appropriate Tow' range including frequen- 
cies lower than 0.5 Hz, and a 'high' range extending beyond 
5 Hz. The exact limits of these ranges vary individually ac- 
cording to AT, T and count rate. The high frequency range 
for XTE Jl 118+480 was not computed as it may be affected 
by telescope vibrations due to wind (Fig. [Tl i Spruit & KanbachI 
2002). For SWIFT J1753.5-0129, on the oflier hand, flie to- 
tal lightcurve duration is too short to produce a meaningful fit 
at low frequencies, while the PSD variance is too low to give 
enough signal at high frequencies, restricting the analysis to the 
base frequency range. 

4. RESULTS 

7. The first important result of Fig.|2]is that XRBs do show 
optical rms increasing with flux over large ranges in frequency. 
Fitting with the linear parametrization of UOl, a = k{F—C), 
the intercept along the flux axis is generally consistent with a 
positive offset. These results are listed in Table |2] 

2. Compared to the X-ray a-F relation, the optical relation 
apparently exhibits more deviations around a linear fit, espe- 
cially towards the highest and lowest flux bins. This is partic- 
ularly evident in the case of XTE Jl 1 18+480 (Fig. |2ji) where 
the scatter in the lowest count-rate bins might even suggest dif- 
ferent behavior in this flux regime. Quantifying the linearity of 
the data with Kendall's rank correlation coefficient (r, as also 
used by G04), the optical shows a weaker correlation (tq) as 
opposed to the corresponding X-ray (tx) values; these are also 
listed in Table |2] 

3. GX 339-4 shows the tightest linear correlation in the op- 
tical with the highest tq value over all frequency ranges inves- 
tigated (also evident from Fig. |2j)). It is also the source with 
the highest fractional rms variability amplitude, as can clearly 
be seen by the fact that its optical PSD has higher broad-band 
noise than the other sources to well beyond 1 Hz (Fig.[TJ. 

4. The rms-flux relation generally flattens (i.e. slope k de- 
creases) with frequency. This trend is consistent with the in- 
terpretation that k changes in accordance with the rms of the 
varying lightcurve component over the frequency range of in- 
terest (UOl). This can be roughly discerned from the PSDs of 
Fig.[T] These have been plotted in units of Frequency x Power 
spectral density, so a horizontal line represents equal variabil- 
ity amplitude per decade. Since higher frequencies typically 
have a lower rms amplitude, the corresponding k values are 
smaller. Furthermore, the lower optical PSD variance (as com- 
pared to X-rays) is manifested in the fact that ko < kx (Ta- 
ble |2]i. Finally, the only case in which the relation steepens 
with frequency is GX 339^, for which fcj^'" < k^^^'^; again, 
this is consistent with the steepening of the optical PSD seen in 
this case. 

5. DISCUSSION 

A linear-Uke rms-flux relation is exhibited by XRBs in the 
optical. The scatter in this relation is apparently larger than in 
X-rays, but this may well be the result of observational biases. 

For instance, the optical data for XTE Jl 118+480 were 
observed without any filters over the full optical wavelength 
range (Table [T]!. The origin of the variable optical emission 
from this source remains unclear, with several components in- 
cluding the disk, accretio n flow and jet l i kely to contribute 
(e.g. iM erloni et al. 200 (1 lEsin et al.l l200lL IChatv et allllool 
lYuan etal.ii2005t iReis et al.ll2009l) . This could allow for the 
emission from distinct components to overlap and introduce 
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Fig. 1. — Optical (red) and X-ray (blue) power spectral densities (PSDs) of each source. The PSDs are normalized so that their integral gives the fractional 
variance of the lightcurves over any frequency interval, and are plotted in units of frequency X Power/Hz so that a horizontal line denotes equal power per unit 
frequ ency-decade. The naiTow spikes above ~ 7 Hz in the optical PSD of XTE Jl 11 8-1-480 may be a re sult of telescope vib rations in high wind tSpruit & KanbachI 
120021) . The optical pealfs below OA Hz for the remaining two sources are real (cf. Gandhi et al. in prep.. [Durant et aU2009t) . 



Table 1 
Optical observations log 



Source 


UT Date 


Telescope/Instrument 


A A (filter) 
A 


AT 
ms 


T 
min 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


XTE ill 18+480 


2000.07.05 


Sldnaltas 1.3m/0PTIMA 


4500-9500 (-) 
5500-7000 (r') 


3 


40 


79.8/63 


2587.5/64 


GX 339^ 


2007.06.18 


vlt/ultracam 


50 


47 


76.8/38 


2712.9/39 


SWIFrjI753.5-0129 


2007.06.18 


vlt/ultracam 


5500-7000 (r') 


39 


31 


28.7/20 


101.5/21 



Col. (4) lists the approximate wavelength coverage (and filter name in brackets); In Col. (5), AT refers to the best time resolution available in the optical, in milli-seconds; Col. (6) lists the 
approximate full length of the lightcurves used in the analysis, in minutes. Cols. (7) and (8) report the chi-square / degrees -of-freedom for the lognormal ('LN') and simple normal ('N') 

histogram fits to the distribution of lightcurve fluxes (see section[5land Fig. [3). 



Table 2 

rms-flux relation fits to the relation (j-k{f-c) 



Source 










Optical 










(I) 


(2) 


(3) 


(4) 


.,0.5-5 

(5) 


(6) 


0.5-5 

(7) 


7.,high 

(8) 


^high 

(9) 


high 

(10) 


XTE Jl 118+480 


0. 11 ±0.05 


1429±11.8 


0.60 


0.09±0.01 


ll4l±2.5 


0.67 








GX 339^ 


0.09±0.04 


-0.82±11.3 


0.87 


0.23±0.01 


12.86±I.I 


0.95 


0.12±0.01 


16.52±0.9 


0.94 


SWTJ1753.5-0127 








0.08±0.01 


18.59±3.16 


0.82 








Source 










X-ray 
















r'r 


j.O.5-5 


^0.5-5 


0.5-5 




^high 


high 


XTE Jl 118+480 


0.23±0.04 


115.99±77.8 


0.87 


O.18±0.02 


llji.5±47.0 


1 








GX 339^ 


0.39±0.07 


10.32±8.8 


1 


0.35±0.09 


17.64±13.2 


1 








SWT J1753. 5-0127 








0.21±0.01 


41.67±17.2 


0.96 









The slope {k) and intercept (C) linear fit parameters are stated for three frequency ranges. The first is a 'low' range including frequencies below 0.5 Hz (Cols. 2 and 3) along with Kendall's 
rank conelation coefficient (r) as a measure of linearity between a and F (Col. 4). This is followed by the conesponding values for the range of 0.5-5 Hz (Cols. 5, 6 and 7) and a 'high' 
range at greater frequencies (Cols. 8, 9 and 10). See Fig.[2]for conesponding plots and details. The fits to the optical are in the first three rows and are denoted by subscript O. The 
coiTesponding parameters for the X-ray fits are below these, with subscript X. The units of C aie cts s~^, with the optical values (Co) being normalized to 1000 cts s~^ for better 

legibility. Errors quoted are for Ax =1 on each fit parameter. 



4 



P. Gandhi 



1.5 



1.0 



0.5 



0.0 



XTEJ1118U80 
Optical 

0.1-0.5 Hz 

0.5-5 Hz 



(a) 



/ A 



5 10 15 20 

Count rate (10^ s"') 



25 



100 
80 
„ 60 
"=> 40 

20 









XTEJ1118+480 
X-ray 



0.1-0.5 Hz 
0.5-5 Hz 



(d) 



100 



200 300 
Count rate (s"') 



400 



500 



5 
4 

3t- 
2 r 

1 i- 
oL 



GX 339-4 
Optical 



(b) 



0.0625-0.625 Hz 
0.5-5 Hz 
5-10 Hz 



♦V 



10 20 30 

Count rate (10^ s"') 



40 



25 
20 
„ 15 
10 

5 






1 1 1 1 1 1 1 

GX 339-4 


' ' ' ' 1 1 




X-ray 






0.0625-0.625 Hz 






- - - 0.5-5 Hz 












— ( 






(e) 






1 ■ ■ ■'^ 1 . . . . 





10 20 30 40 50 
Count rate (s"') 



60 



1.5 



1.0 - 



0.5 



0.0 



SWIFT J1753.5-0127 
Optical 

0.5-5 Hz 



(c) 



10 20 30 

Count rate (10^ s"') 



40 



80 



60 



^ 40 
D 

20 



SWIFT J1753.5-0127 
X-ray 

0.5-5 Hz 



(f) 



100 200 300 400 

Count rate (s'^) 



Fig. 2. — The mis((7)-flux(i^) relation in the optical (plots on left) and X-ray (right). For each source, a common base frequency range of 0.5-5 Hz is plotted 
(green dashed linear fits). Depending on signahnoise, a 'low' range (red dotted) and a 'high' range (blue dot-dashed) were also computed. The exact ranges are 
labelled and the corresponding linear fit parameters are listed in Table|2] 



Optical rms-flux relation in X-ray binaries 



5 



significant scatter in any underlying strong correlation. On the 
other hand, it is to be noted that the X-ray data have lower sig- 
nahnoise than the optical. So the larger error bars could well be 
hiding more scatter in X-rays, because the r statistic does not 
account for measurement uncertainties. 

Furthermore, optical timing observations are sensitive to pre- 
vailing weather conditions. In the case of XTE Jl 1 18H-480, at 
least, there is some evidence of th e wind contaminating the vari- 
ability power (S pruit & Kanbachll2002h . though only the night 
least affected by this bias has been chosen here. Finally, the fact 
that the overall source variability rms levels are much lower in 
the optical is an important limitation. The impact of this is that 
the X-ray relations in Fig. |2] cover a wider dynamic range on 
the X-axis (F) and can hence be better determined than the cor- 
responding optical ones. 

With these caveats in mind, the detection of a linear rms- 
flux relation is important because it provides new constraints 
for any model attempting to explain the underlying mechanism 
of the rapid variability. Uttley and collaborators have ruled out 
additive shot and self-organized criticality models, as these are 
unable to easily reproduce a linear <j-F relation in X-rays over 
a wide range of timescales. From observations of the accreting 
neutron star SAX J1808. 4-3658 and correlating the sour ce rms 
with th e strength of high-frequency coherent pulsations. lUttlevI 
(12004 suggests that the rms-flux relation must originate within 
the accretion flow itself, and not in independent coronal mag- 
netic flares. So-called 'propagati ng perturbation' models are in- 
stead favored (Lvubarskii 1997: lMisralF2000t iKing et al.ll2004t 
iTitarchuk et al.ii2007i) in which accretion rate variations at any 
given radius in the accretion flow propagate towards smaller 
radii to modulate the radiative energy release in all inner re- 
gions. This naturally couples the variability over a wide range 
of radii (hence, timescales) and results in a flux-dependent rms 
amplitude. U05 have shown that this also implies that X-ray 
lightcurves must be intrinsically non-linear; specifically, the X- 
ray flux distribution ought to be have a lognormal distribution. 

The present work now extends these results to the variable 
optical flux of XRBs. We can then also test the above hy- 
pothesis of non-linearity by plotting the distribution of optical 
lightcurve fluxes. These are shown in Fig.[3]for all three XRBs. 
In each case, we fit a symmetric gaussian as well as a three- 
parameter lognormal model (cf. Eq. 3 of U05). The number of 
flux measurements per bin {N) was required to be at least 50 
for the fit, with the bin error assigned as \/lV. 

GX 339^ clearly exhibits a lognormal distribution with a 
very significant tail and even some apparent excesses above the 
fit in the highest flux bins. The other sources are also better- 
described as lognormals, though the difference with respect to 
a simple gaussian fit is more modest, at least to the eye. The 
best-fit values for all fits are listed in Table [T] For complete- 
ness, it is noted that the lognormal count-rate offset parameter 
(r in Eq. 3 and section 4. 1 of U05) was found to be non-zero in 
each case; again, this could be best-constrained in the case of 
GX 339^, where it matches the range of values for the rms- 
flux relation offset (Co). 

As mentioned by U05, lognormal fits may not be formally 
acceptable due to weak non-stationarity of the lightcurves and 
other effects; this may explain some of the excess in the val- 
ues. But, as can be inferred from Fig. |3] the residuals for 
GX 339^ are largely dominated by the highest three (afore- 
mentioned) flux bins alone. So GX 339^ is the source with the 
highest variability (Fig. [U, and it displays the tightest rms-flux 



relation as well as the 'most' non-linearity. Evidence for non- 
linearity is present in the other two objects, but is less striking; 
both exhibit lower broad-band noise. This behavior is exactly 
as predicted by U05, but now applied to the optical lightcurves. 

Whatever the physical mechanism of variable optical emis- 
sion may be (synchrotron from the base of a jet is one possi- 
ble source often discussed in the literature), it should be the 
result of modulations from an ensemble of coupled perturba- 
tions. As compared to auto- and cross-correlation function 
analyses, which are most sensitive to coherent variations with 
the shortest delays, the rms-flux relation reveals underlying 
connections over much wider timescales. The new results pre- 
sented herein are based on only the handful of relatively-short 
optical observations which happened to capture the sources at 
differing stages of their X-ray low/hard states. For instance, 
GX 339^ was observed at the onset of this state with a rela- 
tively faint optical counterpart, while XTE Jl 1 18H-480 was seen 
in its optically-brightest state. The relative fluxes of the vari- 
ous accretion components are bound to change as the sources 
evolve through their states (e.g. iGandhi et al.l2008l and also the 
measurement of the X-ray rms-flux relation in several states by 
G04). Longer lightcurves at multiple optical and X-ray wave- 
lengths, and in various source states, will greatly aid more pre- 
cise comparisons. 

In the low/hard state, the disk is not the sole contributor to 
the optical and X-ray fluxes. Hence the coupling of perturba- 
tions must span not only many radii, but also various accretion 
structures, in order for it to be observed in the C omptonized 
components and the jet (e.g. 'Done et al."2007', and refer- 
ences therein). On the other hand, Zhang (2007) stresses the 
fact that solar flares follow a similar rms-flux relation which, in 
turn, disfavors propagating perturbation models, because there 
is no accretion flow in the solar corona. IZhand suggests that 
a large-scale inter-connected magnetic field in which fluctua- 
tions are able to rapidly diffuse to all parts of the corona may 
work instead. More investigation is needed to understand such 
a model in d etail; but in this con text, we note that the reser- 
voir model of iMalzac et al.l (|2004|) could prove to be consistent 
with all the observations. Their study proposes that the strength 
of the emergent variable power scales in proportion to the total 
amount of energy stored in some 'reservoir'. Though energy 
injection in their model is parametrized as additive linear shots, 
a quasi-linear rms-flux relation is predicted in X-rays; this is a 
consequence of coupling between the fluctuating energy dissi- 
pation rates within the jet and corona. A simple modification 
of the model to include non-linear shot amplitudes (by taking 
the exponential of the lightcurves; cf. U05) will result in a sim- 
ilar relation in the optical as well. This could be achieved with 
minimal changes to the resultant optical/X-ray coherence and 
cross-correlations functions. A good candidate for the reservoir 
is a large-scale magnetic field. Cumulative buildup and dis- 
sipation of energy within the coronal 4-jet (azimuthal-npoloidal) 
field should result in the coupled modulations required by the 
rms-flux relation. 
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Fig. 3. — Binned frequency distribution of mean optical source flux. In each case, a lognormal distribution fit is shown by the red curve, while the dotted green 
curve is the best symmetric gaussian fit. All cases are statistically better-described by lognormals. The y-axes are plotted with a logarithmic scaUng in order to 
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